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1. Introduction 

The main goal of the Hadron Spectrum Collaboration is the calculation of low-lying excited 
hadron spectra. Exploratory calculations of both the Nucleon ^ and Delta spectra have been 
completed as well as preliminary results on two-flavor lattices [|3j] . These spectra may contain 
resonances as well as multi-particle states, both of which are distorted near decay thresholds. All- 
to-all propagators must also be used to create operators that interpolate multi-particle states, which 
consist of single particle states with finite momentum. To extract a large number of both resonant 
and multi-particle states, a large variational basis of spatially extended operators must be used. The 
stochastic construction for all-to-all propagators results in significant computational savings due to 
a complete factorization of the source and sink information. The variance in stochastic estimates of 
all-to-all propagators may be reduced significantly by the dilution method which is employed 
here. While this method has been demonstrated for mesons [^] and simple multiparticle states [Q], 
its effectiveness has not been established for the spatially extended baryon operators required by 
this project. This work is part of an ongoing effort [||] by the Hadron Spectrum Collaboration to 
assess the utility of diluted all-to-all propagators in the extraction of low-lying baryon resonances. 

To this end, the diagonal correlators of three representative baryon operators are calculated 
using various dilution schemes with a single set of stochastic sources. The relative error on these 
correlators is then plotted as a function of the required number of Dirac matrix inversions. By 
examining the expected error falloff from including additional noise sources, it is concluded that 
an increased level of dilution is preferable to multiple time-diluted noise sources. However, a com- 
bination of dilution and an increase in stochastic noise sources may be the most efficient scheme. 
It is also observed that the fractional error for time + spin + color dilution is comparable to the 
point-to-all method. Furthermore, the fractional error on the exact all-to-all is comparable to that 
of time + spin + color + spatial even-odd dilution. 

2. Methods 

2.1 Stochastic Estimation 

The stochastic estimation of the quark propagator, M^^^x^lxo^o) proceeds as follows: 

First, N r random sources, {r]ai(x,f)}, are generated according to some probability distribution. 
The Z4 = {1,-1,/,—/} distribution was used in this work, although similar results were obtained 
using both the Z2 = {1,-1} and U{\) = {e 10 \ 6 £ (0,27r]} distributions. In the distributions 
described above, all elements are given equal probability. After the N r random sources have been 
generated, the linear system 

M {aam (x,t\x',t') <^(x'/) = TlaaM 

is solved for the N r solution vectors, {0« r J(x,?)}. The quark propagator, M^^^x^lxo^o), is 
given in terms of these source and solution vectors: 

M (aa|j3&)( x >*l x °>'o) =£[<M X >0 ^Oo^o)] 

«^L^2(x,0^xo,'o). (2.1) 
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The above quantity is an unbiased estimator of the quark propagator. However, in order to form 
unbiased estimators for hadronic observables, independent random noise sources must be used for 
each quark. For baryons in particular, 3N r independent random sources must be created and used 
to form 3N r solution vectors. 

2.2 Dilution 

The dilution method amounts to partitioning each of the random sources generated above. For 
computational simplicity, a single set of noise sources was used in this work, although the dilution 
method can be used with any number of noise sources. Therefore rather than the N r random noise 
sources described above, a single noise source, T7^(xo,?o)> is partitioned according to a complete 
set of Nd orthogonal projectors, {P^ a ^ b ^(x,t\x' ,t')}. These projectors are applied to form diluted 
noise sources 

N d 

d=\ 

N d 

The Nd diluted solutions are analogously obtained by solution of the following equation: 

M {aalph) (x lt \x'/) <^J(x'/) = 

While any set of projectors satisfying the above properties may be used, simple 'mask' -type 
projectors were employed for this work which correspond to various coverings of the lattice. For 
example the 'time' dilution scheme corresponds to a set of projectors which each have support on a 
single timeslice only. This scheme can be partitioned further to form 'time + spin' or 'time + color' 
schemes, in which each projector has support on a single timeslice and spin or color, respectively. 
Spatial coverings of the lattice may also be used. Projectors in the 'time + spatial even-odd' scheme 
have support on a single timeslice and either the 'odd' or 'even' sites of the lattice, according to a 
red-black checkerboarding scheme. 

The dilution method may also be used in combination with exact solution of low-lying modes 
of the quark propagator [^j. In this case the diluted solutions are projected into the complement 
of the space spanned by the exactly solved low-lying modes. The solution of low-lying modes is 
not used in this work. Furthermore, the use of a single time dilution projector reduces the required 
number of diluted sources to those with support on a single timeslice only. This also reduces the 
required number of Dirac matrix inversions and amounts to forming correlators without averaging 
over the source time. Additional time sources may be added for increased statistics. 

2.3 Implementation for Baryons 

The above method can be used to estimate baryonic two-point functions. This is done by 
forming the 'source' and 'sink' functions, T*f AdBdc \t) and Clf AdBdc \t), each composed of three rj 
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single- singly- triply- 

site displaced displaced 

Figure 1: The three types of elemental extended baryon operators used for this work. Solid circles represent 
smeared quark fields, lines represent smeared link fields, and hollow circles the location of the reference site. 



or (f) fields. These functions are given by 

4d A d B d c ] M _ (£) v c _ r(A)[d A ] , v r(B)[d B 



rr*no =c% Wk Le ab c c w w *%rM x (2.2) 

;(c)[ 



where the c L - .^ coefficients represent both spin and covariant displacement structure for the £th 
operator and the A,B,C indices represent three independent noise sources. Although the operators 
described above have no spatial momentum, this method can be trivially extended to operators 
with finite momentum by the addition of Fourier weights in the spatial sum. These functions are 
then combined via summation over the dilution indices ([dAdgdc]) to form baryonic correlation 
functions. 

In order to access both radial and orbital excitations of baryons, operators which have an 
extended spatial structure must be used. This is achieved by the use of covariantly displaced oper- 
ators in several geometries. These elemental building blocks must be combined to form operators 
which transform according to the lattice symmetries [^J. To assess the effectiveness of various 
dilution schemes, a representative set of three operators was chosen consisting of a single-site, 
singly-displaced, and triply-displaced elemental operator (Fig. [[]). Contamination from higher- 
momentum modes is significantly reduced by using operators which contain smeared quark fields. 
For this work, a Gaussian smearing scheme was employed [jl|] with a smearing radius of a = 3.0. 
The exponential smearing weights were approximated using n a = 32 iterations. Stout-link smear- 
ing [j|] has been shown [jl|] to reduce noise in spatially extended operators. It was employed here 
with n p p = 2.5 and n p = 16. 



3. Results 



To test the above formalism for the baryonic operators described above, 100 quenched gauge 
configurations were used with the following parameters: L s = 12, L t = 48, a s as .1 fm, jS = 6.1, 
and m % pa 700 MeV. The relative error of the diagonal correlator at t = 5 was used as a basis for 
comparison while tests performed using the diagonal correlator at t = 10, 15 showed similar results. 
This quantity is plotted versus the total number of Dirac matrix inversions required (See Figs. ||, 
|3[ and The horizontal axis in these figures represents 1 / \/Ni nv , where Nj m , is the total number 



4 



Stochastic All-to-All Propagators for Baryon Correlators 



John Bulava 



| Single-Site, t = 5 | 









o 

i— 





45 


UJ 






CD 





40 


> 






— ' 

CO 





35 



m 



time 

■ time+space_eo 

A time+color 

▼ time+spin 

# time+color+space_eo 

■ time+spin+space_eo 
^ time+spin+color 

▼ time+spin+color+space_eo 

• pt-to-all 
)K exact all-to-all 

o oo'^ — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — L 

0.00 0.02 0.04 0.06 0.08 

1/V(N_inv) 

Figure 2: The relative error of the diagonal correlator for a single-site elemental operator (at t = 5) versus the 
total number of Dirac matrix inversions required for various dilution schemes with a single stochastic source. 
Along with these dilution schemes, the point-to-all result and the exact all-to-all result are plotted. The dotted 
line represents the error falloff expected from increasing the number of time-diluted noise sources. The solid 
line represents the same for time+color-diluted noise sources. 



rr 0.30 




of required Dirac matrix inversions. In these coordinates, the error falloff due to an increase in the 
number of noise sources for a fixed dilution scheme will be linear with a flattened approach to the 
gauge noise limit. 

As all the subsequent points in these figures lie well below the dotted line, it is clear that 
increasing the level of dilution past time dilution is more efficient than adding time-diluted noise 
sources. However, the improvement is less apparent for the solid line. Here there may be a small 
gain in the fractional error, but computational cost and storage must be considered. An increase 
in the number of dilution projectors by a factor X results in an X 3 increase in the number of 
components in the source and sink functions of Eq. 12. An increase in the number of diluted 
noise sources however, results in the same increase to the number of components of the source and 
sink functions. 



4. Conclusions 

The dilution method for variance reduction has been tested for baryons. As is seen in Figs. ||, 
|3j, and |], increasing the number of dilution projectors usually results in a decrease in error that is 
greater than the naive 1 / y/N expected from increasing the number of diluted noise sources. It is 
also apparent from these figures that time + spin + color dilution possesses comparable error to the 
point-to-all method and that time + spin + color + spatial even-odd dilution is comparable to the 
exact all-to-all result. 
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| Singly-Displaced, t = 5 
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Figure 3: The same as Fig. 2 for a singly-displaced elemental operator. 



| Triply-Displaced, t = 5 
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Figure 4: The same as Fig. 2 for a triply-displaced elemental operator. The black point is slightly offset 
from the pink triangle for clarity. 
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While it is clear that increasing the level of dilution past time dilution is more efficient than 
adding more time-diluted sources, it may not be most efficient to employ dilution schemes finer 
than time + color dilution but rather increase the number of time + color-diluted noise sources. To 
complete these tests, simple fits to the diagonal correlators must be compared, as well as spatially 
extended meson operators. A new method [Q] is also being explored which will yield exact all- 
to-all propagators for a comparable cost. This method is based on a spectral decomposition of the 
quark smearing operator and will retain the desirable property of the factorization of the source and 
sink functions. 
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